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In theories of axion dark matter with large axion decay constant, temperature variations in the 
CMB are extremely sensitive to perturbations in the initial axion field, allowing one to place a lower 
bound on the total amount of inflation. The most stringent bound comes from axion strings, which 
for axion decay constant fa = 10^^ GeV would currently be observable at a distance of 6 x 10^^ 
light-years, nearly 10^ times as far away as our horizon. 
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I. INTRODUCTION 



The Peccei-Quinn (PQ) solution [J, Q to the strong 
CP problem posits a U{1) symmetry which is exact 
up to a QCD anomaly and nonlinearly realized below 
a temperature Tpg. This implies the existence of a 
light pseudo Goldstone boson — the axion [1, 0] — 
whose low energy properties are largely determined in 
terms of its decay constant fa — Tpq and an integer 
Af characterizing the strength of the PQ symmetry's 
QCD anomaly p, H, 0) Hi- The decay constant fa is 
defined so that 9{x) = a{x)/fa is an angle taking val- 
ues 9 G [0, 27r), where a{x) is the canonically normalized 
axion field. To leading order in chiral perturbation the- 



ory, the axion mass is nia ~ 6 x 10~ eV x /-^g , where 
/i6 = (/q/7V)/(10^^ GeV). Laboratory and astrophys- 
ical considerations give /ig ^ 10^^, a bound involving 
model-dependent axion couplings to photons and elec- 
trons 

Cold relic axions from the big bang are a viable candi 
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date for the cosmic dark matter Tl], 11 
symmetry ensures that for temperature Aqcd T < 
TpQ the free energy of the universe is independent of the 
angle 0; as a result, the universe will be populated by 
topological defects in the form of axion cosmic strings, 
around which 9 varies from zero to 27r. Far from these 
strings, 9{x) will assume some random initial value 9i, 
with spatial fluctuations on scales less than the horizon 
damping out with the cosmic expansion. When the tem- 
perature falls to T ~ Aqcd, the QCD anomaly lifts the 
vacuum degeneracy in 9 and there arise M local minima, 
whose zero temperature curvature is given by the axion 
mass. In general, the initial value 9i will not be aligned 
with a minimum of the potential (taken to lie at 9 — 0), 
and coherent axion oscillations will result. For small 9i 
the relic energy density today stored in these coherent 



oscillations is [la, llf 



na ~ n,{^^o^f x 



'5xiovy' 



'16 



3 X lOVi'/' 



/l6 < 1/10 
/l6 > 10 



(1) 
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where He — 0.23 is the cold dark matter that fits obser- 
vation in the ACDM model [17|. The above expression 
ignores nonlinearities in 9i, but one can safely assume 
that for large 9i, the factor of 9f in the above expressions 
is replaced by a number of 0(1). The uncertainty in eq. 
^ and the absence of a formula for 1/10 < /le < 10 are 
due to the crude understanding of the temperature de- 
pendence of the axion mass which arises from nonpertur- 
bative QCD effects when the temperature is T ~ ^qcd- 
In principle, this ignorance could be remedied in part by 
lattice calculations. Nevertheless, it is apparent that (i) if 
9i ^ 0{1) then fa — 10^^ GeV would allow axions to be 
the dark matter today; (ii) if fa > 10^^ GcV then 6'i < 1 
is required; (iii) for /ig ~ 1, interpolation between the 
two formulas in eq. (P) suggests N9i ~ 2 x 10~^ would 
give the correct dark matter density. 

If the universe does not undergo inflation after the uni- 
verse is in thermal equilibrium at a temperature below 
TpQ ~ fa, then causality arguments imply that 9^ ~ 
0(1), which means that one must have fa ^ 10"'^^ GeV 
or else axions would provide a larger value for ildm 
than observed. The ADMX experiment is currently in 
the process of probing this interesting "axion window" , 
IQi" GeV <fa< 10^2 GeV [18]. 

Neutrino masses, the possible unification of the cou- 
plings of the standard model, and the so-called "model- 
independent" axion in string theory [l^, [l^] all suggest 
one should consider a higher decay constant, such as 
fa ~ 10^^ GeV. By eq. ^ this would imply an initial 9i 
which is very small over a region of the universe that cov- 
ers our horizon today, which is possible if inflation occurs 
and the reheat temperature is below Tpg [2l[. Having 
such a finely tuned initial condition for 9i might be ex- 
plained by the anthropic principle [22, [23, [23j [2a| j the ba- 
sic argument being that all possible values of 9i e [0, 27r) 
occur with equal probability in the universe, but only in 
those patches with small 9i are habitable galaxies pro- 
duced. 

We will call this the "ultralight axion scenario": an 
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axion with fa ^ 10^^ GeV, a period of inflation below 
the PQ scale to ensure a homogeneous value of the ini- 
tial axion misalignment angle 9i over today's observable 
universe, and an anthropic explanation for its apparent 
fine-tuning. 

It is important to obtain evidence for or against the 
ultralight axion. In the absence of methods for direct 
detection of such weakly interacting particles, it is natu- 
ral to seek cosmological signatures. Precisely because 6i 
is so finely tuned, the dark matter density today is very 
sensitive to small spatial fiuctuations in the axion field, 
6di{x). In this Letter we discus three different types of 
ultralight axion fluctuations which could survive into a 
post-inflationary Universe and could be detectable today. 



II. AXION FLUCTUATIONS 



positive definite contribution to fla, and hence subject 
to constraints [H, [2^. The argument is that 0^ is re- 
placed by {Of + cr^) in the relic axion energy relation 
eq. ID), where is the fiuctuation contribution. This 
replacement would signify that inflation-induced fiuctua- 
tions make a positive-definite contribution to the relic ax- 
ion energy density today, regardless of how finely-tuned 
9i was prior to infiation. However this is incorrect — no 
matter how much the average value of (9) within a par- 
ticular volume may wander during inflation, after infia- 
tion all values of theta are still equally likely, and there 
are always regions of the universe where {9} is accept- 
ably small — and the anthropic argument then implies 
that our universe today must be located in such a region. 
Therefore the monopole fiuctuation in the axion field has 
no observable consequence. 



Axion fiuctuations are isocurvature, arising before the 
energy density of the universe depends on 9. The ef- 
fect of isocurvature fluctuations on the CMB spectrum 
is greatest at low £, producing acoustic peaks in the 
CMB of much less power than the quadrupole and out 
of phase with what is observed (see [IBl)- Observation 
agrees well with an initial scale invariant adiabatic spec- 
trum; normalizing such a spectrum to the first acoustic 
peak in the CMB, the predicted quadrupole moment is in 
agreement with the COBE result [1^, and hence places 
the constraint that a flat spectrum of isocurvature per- 
turbations comprise < 0.067 of the total perturbation 
spectrum [l7j . This constraint has been applied to in- 
flationary axion perturbations, and results in an upper 
bound on the inflationary Hubble constant Hj for large 

fa ^M,M,M,M,mM,M,M,M,M,m^ we 

recapitulate below. 

The total fluctuation spectrum was fixed by COBE 
[ii] to satisfy {{(ST^))^ w 1.1 x lO^^T , assuming a 
scale invariant spectrum. COBE was primarily sensitive 
to the low £ modes in the fiuctuation spectrum. 

When the quadrupole constraint is satisfied, the £ > 2 
contributions from axion fiuctuations to the CMB are 
small and provide no new constraints. However, it is in- 
teresting to consider the monopole and dipole contribu- 
tions. Dipole isocurvature effects can lead to observable 
consequences - the tilted universe scenario of Turner [s^ 
where the rest frame of the CMB is not the same as the 
matter rest frame. As a consequence, measurements of 
the peculiar velocity of the earth relative to the CMB will 
not coincide with the velocity relative to distant super- 
novae. The constraints on the isocurvature dipole mo- 
ment are currently an order of magnitude weaker than 
the constraints from the quadrupole moment. Thus the 
dipole could only be observable when it is enhanced rela- 
tive to the quadrupole, such as by the geometric enhance- 
ment (fc/fchor)^ which occurs when superhorizon modes 
with k < fchor dominate (37j . This can arise from both 
thermal or topological fiuctuations. 

Monopole fluctuations have been cited as making a 



A. Scale invariant fluctuations 

The first type of ultralight axion fiuctuation we con- 
sider are those induced during infiation, and we repro- 
duce a bound in the literature. During infiation the axion 
field develops a perturbation spectrum for modes cross- 
ing the causal horizon of the form 



27r2 



Hk 
2^ 



(27r)353(k-k') 



(2) 



where is the value of the Hubble parameter when the 
mode of wavenumber k crosses the horizon. For slow roll 
infiation with a nearly fc-independent Hubble parameter 
Hk — Hi, fiuctuations in the axion field eq. ^ lead 
to t > 1 multipole moments in 9i dominated by modes 
with wavelength of order the horizon, with amplitude 
\59\i « Hi/2nfa. 

Following the calculations of ref. [3], the axion con- 
tribution to the total fiuctuation spectrum is 



T 



12 Q.a 59 



(3) 



which, when combined with the COBE and WMAP 
bounds implies 

^^^^<y^{lA^lO-'){i)Ml) = %.2xlO-' . (4) 

Assuming that ultralight axions constitute all of the 
dark matter, and are not diluted by late epoch entropy 
production, the constraint on 59 implies a bound on Hi 



fU'^ X 5 X 108 GeV, /i6 < 10-1 



max ) J16 



Hi < -fff""' = 



„l/4 
J 16 



X 3 X 10^ GeV, /i6 > 10 



(5) 



(For intermediate values of fa we interpolate between 
these bounds.) Infiation at such low scales is inconsis- 
tent with the simplest single field infiation models [ssj . 
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Observation of tensor modes in the CMB from an infla- 
tionary spectrum of gravitational waves will be possible 
in future experiments [3§| provided Hj is high, about 
10"'^^ GeV, which would rule out the ultralight axion sce- 
nario unless fa ^ Mp] such a high value for h as been 
argued to be inconsistent with string theory [201 l40l . l4ll | . 

B. Scale dependent fluctuations 

Evidence for scale dependent axion fluctuations would 
be interesting, giving us information about inhomo- 
geneities preceding inflation. Inflation was discovered as 
a means to make information about initial conditions in- 
accessible to the contemporary observer, but ultralight 
axions — being flne-tuned by the anthropic principle — 
provide a probe that is more sensitive to initial conditions 
than the observed flatness of the universe or the small- 
ness of 5T/T in the CMB spectrum in generic ACDM 
models. The question is: how completely did inflation 
erase information about prior inhomogeneities? 

Two length scales upon which the axion fluctuation 
spectrum could depend are the size of the patch destined 
to become our horizon today, and the thermal wavelength 
at the beginning of inflation. (Other scales can arise from 
non-thermal primordial axion fluctuations — we will con- 
sider those of topological origin below — or from nontriv- 
ial time dependence of the Hubble scale during inflation, 
which we do not discuss). We take iVmin to be the min- 
imum number of e-foldings required so that our entire 
horizon today was causally connected at the beginning 
of inflation, and N to be the total number of e-foldings 
that actually occur during inflation, with TV > iVmin- 
Then the quantity 

Rx ^ e^-^-" (6) 

has a geometric interpretation as the ratio of the horizon 
size at the start of inflation, to the radius of the patch des- 
tined to become our present day horizon. Fluctuations 
of the size of our horizon have comoving wavenumber 

fchor = RxHi, (7) 

where Hj is the Hubble parameter at the beginning of 
inflation. Equivalently, Rx times our current horizon size 
gives the present size of the local patch of the universe 
which was in causal contact at the start of inflation. 

Assuming that the preinflationary universe was not es- 
pecially flat or uniform, observations can provide lower 
limits on Rx ■ For inflation to explain the observed flat- 
ness of the universe, assuming that the curvature radius 
was 0(1) times the Hubble radius at the beginning of in- 
flation, requires Rx — — > 10, given limits on fl 
from WMAP data • To avoid an unacceptably large 
CMB quadrupole moment arising from generic 0(1) pre- 
inflationary classical field fluctuations requires Rx ^ 200 
[sgI ] . Here we are interested in how ultralight axion mod- 
els lead to constraints on Rx, or the existence of an ul- 
tralight axion might be inferred from observation. 



1. Preinflationary thermal fluctuations 

We next consider observable effects due to preinflation- 
ary thermal fluctuations in the axion field, assuming that 
thermal equilibrium existed at a temperature T Tpq. 
As Goldstone bosons, axion interactions fall out of ther- 
mal equilibrium rapidly below that temperature, after 
which their spectrum follows a redshifted thermal distri- 
bution for each mode until it crosses the horizon during 
infiation. To discuss preinflationary thermal fluctuations, 
it is convenient to define the quantity Rt in terms of the 
temperature Ti at the start of inflation, where 

Ti = RtHi. (8) 

Parametrically we expect Hj ^ Tf /Alp. 

In terms of the comoving wavenumber k and cosmo- 
logical scale factor R{t), the correlation between axion 
modes with Hj < fc < T/ is given by 

(SalSa^,) -lL_(2nfSHk - k'), (9) 

We have defined the scale factor R{t) to equal one at 
the start of inflation. For k > Tj Boltzmann suppres- 
sion exponentially damps the thermal fluctuations, leav- 
ing only the zero-temperature scale invariant fluctuations 
discussed in the previous section. 

After inflation has begun, axion fluctuations cross the 
horizon at time tk where k/R{tk) = H{tk)- In models of 
slow- roll inflation H{tk) — Hi and so the thermal fluctu- 
ations that wc eventually see are given by 

((Sak(5ak') ^ -y-(27r)3,53(k-k'), (slow roll). (10) 

Note that these fluctuations are parametrically enhanced 
by a factor oiTi/k > 1 relative to the inflation- induced 
fluctuations in eq. Again, they are isocurvature 

fluctuations and are most important in the dipole and 
quadrupole moments. 

It is convenient to directly compare the contributions 
of the thermal fluctuations eq. Q to low multipole mo- 
ments of the RMS energy density \6p/p\ — 2\6d/9\ to the 
contributions from the zero temperature fluctuations in 
eq. (HI). For £ — 0,1,2 the ratio of thermal to inflationary 
contributions is 



{Sp/p)e,thcim. I y^Rr/Rx Rx < Rt 
{Sp/p)iAnt. ^\{Rt/RxY Rx>Rt ' ^ ' 

For Rt < Rx the thermal fluctuations have a UV cut- 
off corresponding to a length scale Tf^, which is larger 
than the size of our horizon patch, k^^^, explaining why 
their effect falls off as {Rt/RxY = (Tz/fchor)^- For 
Rt > Rx the size of our horizon patch is effectively 
the UV cutoff for fluctuations contributing to large scale 
anisotropics. Therefore the sole temperature dependence 
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comes from the factor of Tj in amplitude of the fluctu- 
ations eq. ©, which enters as -y/T/ in the RMS energy 
fluctuation, explaining the Rt / Rx behavior. 

Thus there is a bound from the quadrupole moment 
on thermal fluctuations analogous to the bound eq. ^ 
on in flationar y fluctuations, found by substituting Hi 
Hiy/Rx/Rx- By the relation Hj Tf/Mp, the bound 
can be expressed as the inequality 

V H^^^ Y Jjmax J ^ ■' 

If, for example, we take /ig — 1/10 and Hi = Hf^'^^ = 
2 X 10^ GeV, then we get the nontrivial bound from the 
quadrupole moment Rx > 3 x 10^. If instead Hi = 
10^ GeV, then the bound is weaker, Rx > 100. 

2. Topological disorder 

Cosmic axion strings formed at the PQ phase transi- 
tion provide a source of nonthermal scale-dependent ax- 
ion fluctuations. As the universe subsequently expands 
these strings intersect and detach loops which decay into 
axions. Causality dictates that there should be Q > 0(1) 
strings traversing each horizon volume at later epochs. 
These give rise to gradients in the classical axion field, 
with the variation A6 across our horizon patch being 

> 2.2^ = ^ . (13) 

Khor Rx 

giving rise to a dipole moment in the energy density of 
order > 2iTQ/{9iRx)- Higher multipole moments are 
suppressed by Rx^- The inequality arises because we 
could be accidentally close to an axion string. During 
inflation and the subsequent evolution of the universe, 
A9 and the dipole contribution to the axion fluctuation 
spectrum is unchanged. The presence of 9i in the denom- 
inator of this expression gives a strong constraint on Rx 
for large axion decay constant. From eq. ([T]) we arrive 
at the estimate for the dipole asymmetry qi 

27rs/j7,y fe.sxioVir, fw<io-\^.. 
''^^[n-J "ji.sxio^Af, Ae>io. ^''^ 

A dipole asymmetry of order 10~^ is observed in the 
CMB temperature [42| , and is usually assumed to be due 



to the peculiar velocity of the solar system with respect 
to the CMB. An isocurvature dipole moment, however, 
causes matter to flow toward the region of higher dark 
matter density so that the rest frames of matter and 
the CMB do not coincide [s^. Therefore one might ex- 
pect our velocity relative to distant supernovae to differ 
from the peculiar velocity deduced from the CMB dipole, 
which is of order 10"'^. Various techniques are used to 
determine the average velocity of galaxies relative to the 
CMB at cosmological scales [13, M, El, H 113, 13, It 
has been estimated that such an effect could be detected 
down to of order 10% of the observed dipole with fu- 
ture analysis of large scale cosmological surveys [i^ . For 
axion dark matter with /ig ~ 10, a constraint on the 
isocurvature dipole qi at the 10~* level would be sen- 
sitive to the presence of a cosmic string more than 10^ 
horizon lengths away. 



III. CONCLUSIONS 

An ultralight axion is a well motivated dark matter 
candidate, provided inflation takes place at a sufficiently 
low scale. In this paper we have discussed how preinfla- 
tionary thermal and topological fluctuations of the axion 
give rise to isocurvature quadrupole and dipole moments 
of the primordial perturbation spectrum, providing con- 
straints on — and potentially providing indirect evidence 
for — dark matter, inflation, the preinflationary universe, 
a very high physical scale, anthropic selection of the dark 
matter abundance, and the explanation for the small elec- 
tric dipole moment of the neutron! Note that similar 
calculations may be done for any light, weakly coupled 
scalar field, such as a string modulus. 
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